Estimation of the repeatedly-projected reduced density matrix under decoherence 
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Decoherence is believed to deteriorate the ability of a purification scheme that is based on the 
idea of driving a system to a pure state by repeatedly measuring another system in interaction with 
the former and hinder for a pure state to be extracted asymptotically. Nevertheless, we find a way 
out of this difficulty by deriving an analytic expression of the reduced density matrix for a two-qubit 
system immersed in a bath. It is shown that we can still extract a pure state if the environment 
^SS^ brings about only dephasing effects. In addition, for a dissipative environment, there is a possibility 

, of obtaining a dominant pure state when we perform a finite number of measurements. 
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PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Dv 



I. INTRODUCTION 



The notion of measurement constitutes a key element in the quantum theory and there is a long history of debate on 
this subject [![■ It is, however, relatively recent that the following fact has been confirmed in a real laboratory that the 
measurement, which is usually described by a projection operator for simplicity but has to be replaced, in a rigorous 
sense, with a generalized spectrum decomposition H, does affect the dynamics of the quantum system drastically 
and in an essential way. The quantum Zeno effect [3(, the effect caused by frequently repeated measurements and 
resulting in hindrance of the dynamical evolution of the quantum system [H just measured, is one of such typical 
and well-known examples and has been discussed quite intensively after its first experimental observation Q. It has 
, also become clear that the action of measurement can have much more profound effects on the quantum systems 

■ than one naively expects or imagines from such a phrase like "wave function collapse." Indeed, the effects have been 
shown to transfer to other quantum systems, not directly measured, but in interaction with the system measured, 

■ and the action of measurement can be utilized to drive even such quantum systems not directly touched to a pure 
state irrespective of their initial mixed states. This is the essence of the proposal of the repeated-measurement-based 
purification scheme [|| , the applicability and effectiveness of which have been examined in various cases Q ■ 

When one considers a possible implementation of such a theoretical scheme to experiments in the laboratories, 
it is necessary to examine its robustness against various imperfections to be encountered in the experiments, which 
may be expressed as a sort of decoherence from the environment. According to the analysis in Ref. @], the projected 
dynamics of the quantum system suffers from the decoherence effect and it is shown that such an effect can deteriorate 
the ability of the purification scheme to such an extent that no (nontrivial) pure state is able to be extracted in the 
limit of infinite number of projections. While this is somewhat an anticipated and discouraging result, the numerical 
estimations of the purity of the state to be extracted for a qubit system, which is coupled with another qubit that 
is subject to repeated interrogations and is immersed together with the latter in a common bosonic bath, show that 
there are parameter regions where states with high purities can be extracted On this basis we have studied 
the possibility to extract a quantum state with a higher purity under a dissipative environment by this purification 
scheme, even though it can not survive in the limit of infinite number of measurements because of a kind of no-go 
theorem Q . This paper is devoted to explore further such a possibility of extracting a dominant pure state at a finite 
number of measurements under a dissipative environment, as well as to point out its robustness against a dephasing 
environment, at infinite number of measurements. 

It is desirable to obtain an explicit expression of the reduced density matrix for such investigations, though in the 
previous study [H the ability of the purification scheme is discussed exclusively in terms of purity. Since solving ana- 
lytically the projected dynamics with decoherence, that is the evolution described by a master equation supplemented 
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with projections, is in general quite difficult, we confine ourselves to a simplest possible case, that is, two mutually 
interacting qubits, immersed in a common bosonic bath. We analyze the reduced density matrix of one of the two 
qubits when the state of the other qubit is regularly projected to one and the same state, according to the spirit of 
the purification scheme [||. It is shown in the next section that if the bosonic environment causes only dephasing 
effects on the qubits, the purification scheme is robust enough to ensure one to extract a pure state, provided an 
appropriate state is confirmed on the other qubit at suitable intervals, just as in the ideal case. On the other hand, if 
the environment brings about a dissipativc effect on the two-qubit system, it is quite difficult to solve the projected 
reduced dynamics in a compact form. However, in Sec. IIIII it is shown that one can still estimate the behavior of 
the state analytically for a large, but finite number of measurements N. After an explicit expression of the projected 
density matrix for one of the qubits is derived in a closed form in Sec. IIII A[ its asymptotic behavior for large N is 
presented in Sec. IIII B[ Even under a dissipative environment, one can find a possibility of extracting a dominant 
pure state when the number of measurements N is large, but not extremely large, provided the dissipation is weak 
enough. However, it is to be noted that in the limit of infinite number of measurements, only a mixed state, except 
a trivial pure (ground) state, can survive, in accordance with the general theorem shown in [8[. We close the paper 
by giving the summary and discussions in the final Sec. IIVI A few Appendices are added to clarify some details that 
arc not shown in the text. 



II. PURIFICATION UNDER DEPHASING ENVIRONMENT 



Consider a quantum system composed of two mutually interacting two- level systems (qubits), X and S, immersed 
in a common bosonic bath. Assume that the dynamics is given by the total Hamiltonian of the form 



Ht = Hxs + -ffbath + -ffint , 

where Hamiltonians for the two-qubit system and for the bosonic bath read 



^xs = |af) + ^f) +5 (afV S) 



Hbath = J dkLo k a\a k , 



and the interaction Hamiltonian 



H 



int 



: J dk (h{k) 



cik + h.c. 



(1) 

(2) 
(3) 

(4) 



causes dephasing on the system. Here the notations are standard, cr£ = (ci ± ioy 1 )/2 (i = X, S) etc. and £ = 
er( x ) + cr( s \ and the interactions with the bath are assumed to be symmetric between X and S and are given by the 
same function h(k), for simplicity. Since the above Hamiltonian Hxs is diagonalized, 



Hxs = E 2 \2)(2\ + E Q \0)(0\ + E + \+)(+\ + E-\-)(-\ 



where 



|2> = |T)x ® |T)s = ITT), 
|o) = |U), 

l±) = -L (yi±(u>-n)/2E+\n) ± y/it^-^Ae+iit}) , 

and the interaction is rewritten in terms of these eigenstates as 



E 2 = -(w + fi) 



(5) 



(6) 



e+ = \/ z (w-n)*+ 5 2 



H int = (|2)(2| - |0)(0|) J dk (2h{k)a k +/i.c), 
the usual procedure Q yields the master equation for the system 

Pxs = -i[H xs , p xs ] + 7 [(|2)(2| - |0)(0|) P xs(|2)(2| - |0)(0|) - | { 1 2> <2| + |0)(0|, p XS } 



(7) 



(8) 
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with the Hamiltonian Hxs , which is still diagonalizcd as in (J6]) but with an energy shift AE, 

E 2 =E 2 + AE, E Q = E a + AE, E± = E±. (9) 

The (dephasing) constant 7 and the energy shift AE are given by the spectral density as usual @. A remarkable 
point is that the structure of the master equation is not affected by the temperature as in the dissipative case (See 
Sec. IIII[) . and the tem per ature only influences the values of the decay constant and of the energy shift. The master 
equation JHJ) is solved [ldl ] in the following form [ll[ 



Pxs(t)= K lPxs (Q)Kj, 

i=0,± 



(10) 



where the three operators are given by 

K = e -^ 2 *- 7 * /2 |2)(2| +e-^°*-^/ 2 |0)(0| + 



e "H+X+ 

K+ = v/cosh 7 t- 1 ( e - ,; ^ 2 *- 7 * /2 |2)(2| +e-^°*- 7 * /2 |0)(0| 
K_ = v /sinh 7 t (e- l ^ t -^ t/2 \2)(2\ - e ~^o* — ^*/ 2 |0><0| 



-)(- 



(11) 



Observe that these operators are all diagonal and the environment causes no transitions between different levels. Only 
the phases between different levels are affected by this type of interaction with the environment and this is the reason 
why it is called a dephasing interaction. 

According to the spirit of the repeated-measurement-based purification scheme Q , we measure qubit X to confirm 
that it is in the same, particular state, while we do not touch the other qubit S, even though it is affected indirectly 
through its interaction with X. This process is repeated many times at regular intervals and we are interested in what 
happens to the state of qubit S, which also suffers the dephasing effect from the bath in addition to the projective 
actions on X. ft has been shown in Q that for a pure state to be finally extracted in this kind of process with 
projections and decoherence, it has to be one of the simultaneous cigenstates of all the relevant projected operators 
and furthermore its eigenvalue for the projected dynamics has to be the largest in magnitude, preferably close to unity. 
In our case of dephasing environment, it is easy to find such simultaneous eigenstates of the projected operators, as 
will be shown below. 

Let us consider to measure the up state ||)x regularly at t = tit (n = 1,2,.. .),since we already know that this type 
of measurement results in an optimal purification of qubit S in the ideal case [7[. The relevant projected operators 
Vt(i) = x(t|-Ki|T)x read (the subscript s shall be suppressed in the following) 



Vt(o) = x(f|A" |T)x = e 



-iE 2 T--yT /2 



I TXT I 



uj — 

cos E + t + i—r^: — sin_E+T 



2E 4 



14X41, 



^T(+) = x(T|A + |t)x = V cosh 7'' 



le 



-iE-2T — JT j'2 



ITXTI, 



V n _ } = x(T|A_|t)x = v^nh^e^ 2r - 7T/2 |T><T|. 



(12) 



It is clear that both the up and down states (of qubit S) are the simultaneous eigenstates of all the above projected 
operators and the eigenvalues and eigenstates for the projected dynamics, which is given by the following map 



p(kr)= J2 V Ul) p({k-l)T)vl %)1 

i=0,± 



(13) 



are read from 

E 

i=0,± 



ITXTI^) = (| e -^-^/ 2 | 2 + |Vcosh 7 r-le- 



i-E 2 T-7T/2| 2 



+ | v ^mh 7 7e-^^-^/ 2 | 2N ) ITXTI 



[(<■; 



14X41^ 



1(0 



ITXTI, 



cos E + t + i- 



i=0,± 

E 

i=0, ± 



2E A 



■ sin E + t 



14X41 = I£I 2 I4X4I, 



E % 

i=0, ± 



%)ITX4|F/ W 
I4XTI4, 



_ -i-E2T-7r/2 



lu — Q 

cos E + t — i— — — sin E+t 

2Ea- 



u> — O 

cos E + t + i—r-^ — sm E+t 



ITXTI 



-iEo,r— J2£* 



riT)(4i, 



2E 4 



|4)(TI=e 



_ p iE 2 r- 1T /2 



CUXTI- 



(14) 
(15) 

(16) 

(17) 
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Starting from a factorized initial state, pxs(O) = |T)x(t| ® p(0) [T^] ? the density matrix is given by 



p(Nr) 



( 



ptt(°) 



e -iE 2 T- 7 T/2 e \ pu(0) 



N 



,iE 2 T--yT/2, 



N 



(18) 



Kr pu(o) 



Observe that the up state |t)(T| is the eigenstate belonging to the eigenvalue unity and that there are no 7's left 
in its relation (fT4|) . It means that the dephasing has no effect on the ability of the purification scheme. The other 
cigenstatcs belong to eigenvalues smaller than unity if sin E + t 7^ 0, just as in the ideal case with no decoherence. The 
above relations clearly show that we will be able to purify qubit S to the up state with no loss of probability, unless 
cos E + t = ±1, irrespective of its initial mixed state, when the qubit X is repeatedly confirmed to be in the up state. 
The repeated-measurement-based purification scheme is thus shown to be robust enough against the dephasing effect, 
at least when |T)x is measured. On the other hand, it is easy to demonstrate that the only two possible pure states 
that can be extracted under dephasing environment are ||) and ||). Indeed, the projected operators x{4>\K± \4>)x, V<^> 
are always diagonal in the basis {|t), ||)}, and the operator Kg admits one of such two states as an eigenstate only 
when |T)x or ||) x is measured. It is worth stressing the possibility of extracting a pure state, i.e. ||), even under 
dephasing in certain conditions. 



III. PURIFICATION UNDER DISSIPATIVE ENVIRONMENT 

Consider next a dissipative environment and assume that it interacts with the two qubits, X and S, through the 
interaction Hamiltonian 



H int = J dfe/i(fe)(E+a fc + E_4), S±=4 X) +4 S) - 



(19) 



For the sake of simplicity, we shall set uj = f2 [13[, though generalization would be straightforward. The system 
Hamiltonian reads 



i/xs = ^3+. 9 (4 x y s) +4 x) CT f)= £ Ei \i)(i 



(20) 



i=2,0,± 



where eigenenergies and eigenstates are somewhat simplified 

|2> = ITT), 
|0) = |U), 

l±> = ^(lU>±IIT>), 



— — — Eq, 

E± = ±9 



(21) 



and the interaction Hamiltonian is rewritten as 

H int = J dkV2h(k) [(|0)(+| + |+)(2|)4 + /i.c. 
When the bosonic bath is at temperature T = 0, the master equation is derived under the usual conditions 



(22) 



Pxs - -i[H XS , Pxs] + 7 [l+)(2|pxs|2)(+| - £{|2)<2|, p XS }] + T [|0> <+|pxs |+> <0| - ^{|+}(+|, pxs}]- (23) 

The decay constants 72^+ and 7^ > (we assume Q > g), given by the on-shcll form factors, arc assumed to be the 

same 7 = 7 2-^+ = 7h »0: just for simplicity. The Hamiltonian Hxs is still diagonalized by \i) (i = 2,0, ±) in (|2ip 



^xs- E 2 = E 2 + AE 2 , E+ = E + +AE+, E = E , E_=E_. 



(24) 



i=2,0.± 



We follow the recently developed technique [lCj to obtain the solution of the master equation (f2"3"| in the following 
compact form 



pxa(t) = e A Vxs(0)e At * + (l-e^^OSo / 9xs(0)Boy(l-e" 7 *- 7 te- 7 ') J BoBi /9 xs(0)Bi t -Bo t +7te" 7 *S 1/9xs (0)B 1 t , (25) 
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where 



e -<^t-7*/2|2)(2| +e-^+*- 7 */ 2 |+)(+| +e- iEot \0){0\ + 



-iE-t 



l-X-l 



(26) 



and 



So = |0)(+| 



Si = l+><2|. 



(27) 



Notice that one can think of a similar master equation that describes the dynamics of the two qubits X and S, in 
interaction with a common bosonic bath at temperature T = through the same coupling (|19j) . but with a different 
mutual coupling from For our purpose, however, it is desirable to make things as simple as possible, for, as will 
be seen below, even the above seemingly simplified dynamics (|23p can bring us with quite involved expressions for the 
reduced density matrix of S when the state of X is periodically projected on one and the same state. Therefore, in this 
paper, we exclusively consider the dynamics (|23p and endeavor to disclose the asymptotic behavior of the projected 
reduced density matrix for system S. 

If the system X is to be repeatedly measured at t = tit (n = 1,2, . . .) to confirm that it is in the state |a)x = 
Q ; |T)x+ |l)x (the normalization factor (1 + |a| 2 ) -1 / 2 is tentatively omitted here for notational simplicity), the relevant 
operators for the dynamics of qubit S read 

xHe A » x = e^-^/ 2 H 2 |T)(T| + \e-^-^' 2 (a*\[) + |T>) (a(j| + <T| 



+ e-^IDUl + \e~ iE - T {a*\i) |T>) («UI - (Tl) 



| a |2 e -i_E 2 r- 7 r/2 + Ig-iB+T— yr/2 + \ & -iE-T 



(c"^ +T ~ 7r/2 + e~' lE - T ) + e~ iEoT 



V, 



(a|B |«)x = ^|i)(aai + (T| 









1 


o 


72 


V2 



Co, 



(a\B B 1 \a) x = a\l){\\ 




a 



x(a|Bi|a) x = ^(a*|i) + |T))(T| 



n 

LV2 







c 2 



(28) 
(29) 

(30) 
(31) 



The projected reduced density matrix p(r) of qubit S is given by the following map 

p(r) - Vp(0)V^ + (1 - e-T T )C 0/ 9(0)C t + (1 - - 1T e-^)C lP {Q)C^ + 7 r e -^C 2 p(0)C 2 t . 



(32) 



A. Projected reduced density matrix 



It is not difficult to confirm that the down state \i) can be a simultaneous eigenstate of all the above relevant 
operators if a = 0. This corresponds to the case where qubit X is confirmed to be in the down state and qubit S is 
projected to the down state ||)(|| with probability 1. This fact has already been pointed out in Ref. @. Furthermore, 
one can easily see that for nonvanishing a, these relevant operators do not have a common eigenstate and therefore 
no pure state can be extracted in the limit of infinite number of measurements, according to the general theorem [8[. 
Nevertheless we are interested in the asymptotic behavior of the projected reduced density matrix p(Nr) after N 
successive measurements, in the hope of finding a way out of, or bypassing such a "no-go" theorem. 
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In order to find the asymptotic form of the projected dynamics when the projection is repeated many times, observe 
first that the above operators Co, C\, which are expressed as (now the normalization constant is recovered) 



Co 



1 



(1,«), Ci 



a fO 



V2(l + H 2 ) V 

in the up-down basis, enforce the system to be in a pure (down) state 

'0 



1 + lal 2 VI 



(1,0) 



J) (o,i) = U)UI = pu, 



while the operator C2, which can also be expressed as 

c 2 = 

drives the system to another pure state 



\/2(l 



|2-i a 



(1,0), 



I 



1 



1 + \a\ 2 \a 



* (l,a) = I a*) (a* I = p a *, 



(33) 



(34) 



(35) 



(36) 



irrespectively of the state the system had lived in just before projection. This means that in general the projected 
dynamics drives the system to a mixed state, which would be made more apparent if it is written in the following 
form (apart from normalization) 



p(r) = Vp(0)V^ + F(p(0)) Pu + G(p(0))p a ,, 
where the positive constants F and G read 

F iPm = }~ e ~ZA a*W)\<**) + „ '"^ (l - e- JT 7re-^)(T|p(0)|T) 



2(1 + |a|») 



and 



G(p(0)) 



2(1 + |«P) 



7 re-^(T|p(0)|T), 



(37) 



(38) 



(39) 



respectively. 

It is not difficult to see that after N repetitions of the projected dynamics, the state of the system S, apart from 
the normalization, is driven to [po = p{0)] 



N-l 



Pn 



p(Nt) = V N p(0)V^ N + E [FfaW-^PuVW- 1 -" + GipJV^-WVi"- 1 



-i-fci 



(40) 



fc=0 



The coefficient F(pw) satisfies 
1 - e-T T 



F(p N ) 



2(l + l«| 2 ) 

1 - e"T T 

2(l + l«l 2 ) 

iV-l 



a*\p N \a*) + 



|a| 2 (l-e-T T -7Te-^) 



(TIPivIT) 



E F (^) 

fc=0 

E 



(a*\V N p(0)V^ N \a*) + 
1-e-T 



(1 + H 2 ) 2 

|a| 2 (l- e-^ T -7re-T r ) 



L2(l + |a| 2 ) 
1 - e -Tr 



a*V 



N-l-k 



(1 + H 2 ) 2 

2 M 2 (l- er~> T -7re" 7T ) 



fc=0 



2(1 + H 2 ) 



[a*\V 



N-l-ki 



(l + |a| 2 ) 
2 |a| 2 (l-e-T r - 7 re-T r 



m N p(0)V^ N n) 

W-l-fci. 



2\2 



(l+l«| 2 ) 2 



(T|F J 



(41) 
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Similarly, wc have 



G(PN)= l l ^tm (TlpivlT) 



2(1 + |a| 2 ) 
|a| 2 7Te _7r 
2(l + l«| 2 ) 

N-l 



(T|^(0)^ t7V |T> 



JV-l 



a\ jre 



-7T 



(i + H 2 ) 



(Ti^-^ia*) 



(42) 



It is now clear that these coefficients satisfy the recursion relation of the following form 

N-l 



'F(Pn) 
G( PN ) 



E^ 

fe=0 



N-l-k 



G( Pk ] 



>N ■ 



(43) 



where two-by-two matrices An-i-U and two-component column vectors 6^ would be evident from the previous 
expressions. If a parameter x is introduced, its formal solution is easily found through an identity relation 



F(pn) 
G(Pn). 



(44) 



fe=0 



that is, the coefficients F(pk) and G(pk) are given as those of the fcth power x k on the right hand side. Plugging these 
values into (|40[) . wc obtain the explicit expression of the density matrix of qubit S at Nr. 
We assume the diagonalizability of V and introduce the eigenvectors of V 



V\ui) = \i\ui), (vi\V = Xi(vi\, ^2 \ui)(ui\ = 1, (vi\uj) = 5 K 



(45) 



i=l 



to expand it as 



7 = ^AiK)(t 



i=l 



The elements of the two-by-two matrix Ak and the column vector bk are then expressed as 

a, b a,b 



where 



A a fc = A a A b 



(46) 



(47) 



(48) 



and the explicit expressions of (C ab )ij and (d ab )i are found in Appendix [A"l The above identity relation (|4"4")l tells us 
that the coefficients F(pk) and G(pk) are explicitly given by 



'F(Pk) 
G(Pk) 



E 



Q<£,k 1 ,...,kg,n<k 

i + k 1 + ... + k i + n =k 



k k-i k-i-k-L k-t-kx A^„i 

EE E ••• E E (A ai& j fci ---(A a ^)^(A af+A+1 ) fc -^-----^c^^...c^^ 

£=Q k 1= k 2 =0 k e =Q "i,t>j 



(49) 



Observe that the indices k\,...,ki only appear as the exponents of the eigenvalues. The summations over these 
indices are not difficult to perform and each summation results in a similar expression. Actually, since we can prove 
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by induction (see Appendix [B)) . 

k-i k-i-k x k—e—ki—— 

EE- E 



k± k,£ k—l—ki kg 



X i ' * * X q X 



fe 1= fc 2 =0 



k e =0 



J i+1 



i+1 



(xe+i - xi)(x e+ i - x t -i) ■ ■ ■ (xe+i - xi) 



(x e - x e+ i)(x e - Xt-i) ■■■(X£- X\) 
x\ 



(xi - xe+i)(xi — xe)-- - (xi - X2) 



t+i 



E 



^ Y\ {-^m X n ) 



l<n<t + l 



we arrive at 



F(pk) 
G{p k ). 



k l+l 

EE E 

e=o ro=l *u,i>p 



(K mb J k 



i< a ,0<e+i !<„^ <+1 

Notice that pn in (|40|) is expressed in terms of the eigenvectors |ttj) and \vj) as 

2V-1 



F( P k)' 
G(Pk). 



(50) 



(51) 



(52) 



i,j i,j k—0 

Plugging the explicit forms of F{pk) and G(pk) (|5Tj) into (|52|) . we are left with the following terms to be evaluated 

iV-l 



-i-fe 



fc=0 



'F(pk) 
G(pk) 



N-l 



E(A,r-TE E 

k=0 £=0 m=l aafip 



(K mbm ) h 



II ( A a„,6 m - K n b n ) 



Qaibi _ _ ^afbf^af + ibf + i 



(53) 



l<"./3<*+l 1<„< <+ 1 



The summations over k and f can be interchanged ^2^=o Sfco = Sfco 1 EfcL/ ano - that over k is performed to yield 



jv-i e+i 

EE E 

£=0 m=l "o,^ 

1<o,/3<* + 1 



II ( A a m 6 m -A Q „ 6 J 



(Ajj) (A anl b m /A;j) (Aq m fc m ) ^.a-^ ^a f 6f ^a e+1 b e+1 



Lin ~ A„ 



(54) 



This is the explicit form of the second term of the rhs of (|52|) and brings us with the exact expression of pn = p{Nr). 



B. Asymptotic behavior of the projected reduced density matrix 

We are now in a position to evaluate the asymptotic form of for large N . The particular form seen in (|54[) 
implies that the dominant contributions are due, in general, to those terms with the same set of indices {ij} = 
{a\bi} = {0262} = • • • = for in such a case there are order-A^ terms contributing constructively. This 

can be understood by looking at the original expression (|49|) . See Appendix [C] for more explanations. The dominant 
contributions to the second term on the rhs of (j52|) 



N-l 



i,j k=0 



'F(pk) 
G{ P k) 



EE^)^" 1 ""!^)^!^!^!^), (v<\Pa-\Vj) 

= \ul}{ U A(( V l\PU\ V l)' (vi\Pa*\vi)^Mll + \u 2 )(u 2 \(^(v 2 \pu\v 2 ), (v 2 \p a *\v 2 ) ) - 

+ \ui)(u2\((vi\pa\v 2 ), (vi\p a *\v 2 )^M 



h.c. 



(55) 



9 



are found in the above (column) vectors M.if. 

JV-l JV-l JV-l k-l k—e—ki 

Mn~ E<r^o + E EE- E 

«=0 C=0 fc=J>+l fei=0 fc^=0 

JV-l k-l 



N-l-Z /y?ll\£ jH 



(A n - A)~ 1 (Af 1 - 06„ + e E ; as-^cc 11 )^ 1 , 



fc=l c=o 



(56) 



the second term of which is evaluated to be 

JV-l k-l 



EE 

fc=X £=0 



A w-i-^ c iiy d n 



1 + a9 



1-1 



An 



JV-l 



d n - Afr'c 11 i-7— 

An 



1 - 



An 



JV-1 



(57) 



Similarly, wc have 
N-l 



N-l N-l k- 



k-i-kx 



m 22 ~ £ a£- x ^6o + E E E • 

e=o i=o k=i+i fci=o 



E 



A 



JV-1- 



t=0 



fc-£ 



if) 



A 



JV-l 



C 11 A 2 



Au 
A 22 



JV-l 



U \iv-i 



1 + 



An 



- 1 



c_ 

A 2 2 L 



1 - 



9— 

A 22 



AT-l' 



(58) 



and 



Mi-, 



N-l N-l N-l k-l k-t-ki- 

EAS-^fto + E E E 

£=o f=o k=e+i fei=o 

(Ai 2 -A)- 1 (Af 2 -Ofo 
C 12 



Af^C 12 )- 1 



1 



A 



12 



1-1 



9lL 

Ai 2 



E 

fc»=0 



JV-l' 



A 



JV-l-f //7l2s^ j12 



Af 2 " 2 C 12 ( 1 - 



A 



22 



Al2 



JV-l' 



(59) 



For weak dissipation 7T < 1. the matrix elements of C ab (and _4o) ar e of the order of 7T or higher [see (|A1[) ]. while 
the maximum (in magnitude) eigenvalue Ai is expected to be of order unity Ai ~ 0(1) with corrections of order jt. 
In this case, the dominant contributions to Mij are estimated to be 



JV-l .11 



Mn ~ (JV - l)K d 



M 



22 



AffV 1 , A4 12 ~(JV-1)A 



JV-l .12 
12 " ■ 



Therefore, if the above pjv (|52p is suitably normalized, the state is approximated as a pure state 

p N ~\u 1 )(u 1 \ + 0(< J T, |A 2 /Ai| JV ), 



(60) 



(61) 



for a large, but not extremely large N under which the above approximation is valid. 

Notice that the validity of the above expression (|60|) and therefore that of the ensuing relation (f6"Tj) are limited. 
Actually for a larger N, N^t 3> 1, the approximation (1 + x) N ~ 1 + Nx with x = 0(jt), on which the expression 
(f6"0"|) has been based, is no longer valid. Instead, since 



*JV _ e JVlog(l+a;) ^ e N(x-x 2 / 2 +---) 



(62) 



we expect, for a larger iV, say I/7T < N < 1/(jt) 2 , 

Mn-Afr'e^ 1 ^^ 11 )" 1 ^ 1 , M 22 ~ Afr 1 ^" 1 ^ 11 , AI12 ~ A^~ 1 (C 12 )~ 1 e (iV_1) ^'d 12 . (63) 

Even though the coefficient Ai\\ multiplying the term |iti)(iti| gives the major contribution, the others coefficients 
are of the same order as the former and never decrease as we increase N . It is easily understood that a similar 
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expression can be found for an even larger A. That is, we are unable to reach a pure state by increasing the number 
of measurements (projections) N, which is in accord with the theorem shown in Q. 

A few comments are in order at this point. The pure state approximately extracted above in (|61|) for an 

intermediate N is not an eigenstate of the projected dynamics ([37]). Actually it is possible to prove that no pure state 
can be an eigenstate of such a positive dynamical map, provided that it includes more than two (in the case of two 
dimensions) projections that (are necessarily not orthogonal and) project to different pure states. Nevertheless, it is 
still possible to consider such a situation where effects of all but one such projections decrease more quickly than that 
of one particular projection and as a result for an appropriately large A the dynamical map effectively becomes single 
dimensional, i.e., an effective manifestation of purification. This is the case in (|6ip . Second, in the strong-dissipation 
limit (7 — > 00) many simplifications are expected to occur in the general form, but the essential points are intact. 
The density matrix p^ takes exactly the same form as the rhs of (|52p supplemented with (|55|) . It is easily seen that 
G(pk) = 0> Vfe and the column vector appearing in ([55]) . C aibl ■ ■ . C aebe d ae+lbe+1 , has only the upper component 



. . . Q a tb[ d a t+l 



(C aibl )n ■ ■ • (C a ^)ii(^ +lf,£+1 )] 



(64) 



The asymptotic form of the density matrix is easily evaluated 

PN ~ wiAiiVXuii^iipuh), H^.h))^ 11 )^ 11 

+ N\\ 1 \ 2 \u 1 )(u 2 \((v 1 \ Pll \v 2 ), (v^lv^i^Y^d^ + h.c. 

+ (N-l)\X 1 \ 2 \u 2 )(u 2 \((v 2 \p u \v 2 ) 7 (vzlp^lv^^y-'d 11 + (N -l)^ 2 ^, p^iC 1 ^-^ 



MMHPuhHOii (Oi + KKHHmMrTi (Oi + h.c 



U2 )( W2 |(«2iPui^)(c 11 )fr 2 (rf 11 )i + pu(c 11 )fr 2 (rf 11 



(65) 



If we neglect differences between (C n )n and (C 12 )n and (d n )i and (d 12 )i, the above expression tells nothing but 
that the state relaxes to the down state pn = On the other hand, if the matrix elements, (C n )ii and (d n )i, 

are greater than the others, (C 12 )n and (d 12 )i, the formers overwhelm the latters as A increases and the density 
matrix will relax to a mixed state 



p N — ► |ui)(ui|(wi|p u |i;i) + \u 2 )(u 2 \(v 2 \p u \v 2 ) +p n . 
If the opposite case were possible, we would have been given another matrix 

\ui)(u 2 \ + \u 2 ){ui\ 



Pn 



(U 2 \Ul) 



MIR 



(66) 



(67) 



but it turns out that this is impossible, for this matrix can be shown to have a negative eigenvalue (the trace of 
the square of this matrix is shown to exceed l!) and therefore it will never show up in our physical process. In all 
cases with strong dissipation, no nontrivial pure state is shown to be extracted, even at an intermediate stage with 
a finite number of measurements. Notice that these observations are actually consistent with the exact result; in 
the strong-dissipation limit 7 — > 00, the operator C 2 becomes irrelevant [see (|32p ] and since the operators Co and 
Cl project qubit S to the down state, the simultaneous eigenstate should be the down state (a = 0), for which the 
operator V becomes diagonal. This means that |«i(2)) - * li(T)) an d Pn — > Pn- 

Finally, it is to be noticed that the dissipative environment also affects the yield, that is, the probability of obtaining 
the target pure state. Actually, the probability of obtaining the above pure state (|6"l"T) is considered to be almost equal 
to the success probability, which is given by the normalization factor of pn in (|52p . if the latter is dominated by the 
pure state as in (|6"l"T) . and this factor is of order ~ (1 + A^7r)A^ for an intermediate A, e.g., 1 <C A < I/7T. On the 
contrary, the dephasing environment causes no effect on the yield. 



IV. SUMMARY AND DISCUSSIONS 



We have shown in this paper that the repeated-measurement-based purification scheme is robust against a dephasing 
environment and the up state of qubit S, |T)(TL can be realized asymptotically without any loss of probability, when 
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qubit X in (rotating- wave) interaction with the former is repeatedly confirmed to be in the up state, |T)x, just as in 
the ideal case. In order to reach the final pure state with fewer steps, we just adjust parameters, say, the interval 
between measurements r, and minimize | cos_E + t + i{ui — fl)(2E + )~ 1 sin E+t\, 

a (lu-n) 2 2 4o 2 2 (w - ft) 2 

cos 2 £ + t + v ' sin 2 E+t = 1 - — - sin 2 E+t > . ^' A7 , (68) 

+ 4El + (u- fl) 2 + 4g 2 + ~ (cj - ft) 2 + 4g 2 v ; 

the equality of which is attained when cosi?+r = 0, realizing an optimal case. This condition is the same as that for 
the ideal case. 

The reason why the dephasing does not affect the ability of this kind of purification scheme may be understood in 
the following way. The dephasing environment surely disturbs the phases of both qubits but causes no transitions 
between up and down states. On the other hand, our purification scheme is dependent on the probability of finding 
a quantum system in some definite state and the change in its phase has no relevance to this scheme. This is a naive 
interpretation of why the ability of this purification scheme is not affected by the presence of dephasing environment, 
when the state |f) is measured. 

We have next considered the case of dissipative environment and examined the ability of the purification scheme, 
though it is already known Q that no (nontrivial) pure state can be extracted in this case according to the general 
theorem. When the qubit X is repeatedly confirmed to be in a definite state, the target system S, in (rotating- wave) 
interaction with X, is forced to be in a definite state, while the surrounding environment constantly drives the system 
to the equilibrium. It seems that a kind of competition between two tendencies, one forced by the projection and the 
other relaxing to an equilibrium state, results in an approximate extraction of a dominant pure state at an intermediate 
stage with a large but not extremely large number of measurements. For a weak damping case, we are able to find 
an asymptotic expression of the state of qubit S, which shows that the dominant contribution is given by one of 
the eigenstates of the operator V (belonging to the largest eigenvalue in magnitude), which is a kind of projected 
evolution operator incorporating partly dissipative dynamics and reduces to the usual unitary operator supplemented 
with projection when there is no dissipation. Of course, if we repeat the measurement indefinitely, the dissipative 
dynamics would swiftly overwhelm the effect of projection (purification) and the system would never be driven to a 
(nontrivial) pure state, which is in accord with the general theorem. 

It is worth mentioning that in actual situations with various causes of decoherence, one has to be careful about 
and aware of the presence of such nonidcal elements and should not perform measurements indefinitely, for there is 
a possibility that an optimal result can be attained at a finite number of measurements, as in our simplified model, 
even though a general criterion seems to be quite difficult to be obtained at present. 



V. ACKNOWLEDGMENTS 



This work is partly supported by the bilateral Italian-Japanese Projects II04C1AF4E on "Quantum Information, 
Computation and Communication" of the Italian Ministry of Instruction, University and Research, and 15C1 on 
"Quantum Information and Computation" of the Italian Ministry for Foreign Affairs, by a Grant for The 21st Century 
COE Program (Physics of Self-Organization Systems) at Waseda University from the Ministry of Education, Culture, 
Sports, Science and Technology, Japan, and by a Grant-in-Aid for Scientific Research (C) (No. 18540292) from the 
Japan Society for the Promotion of Science. Moreover, the authors acknowledge partial support from University of 
Palermo in the context of the bilateral agreement between University of Palermo and Waseda University, dated May 
10, 2004. 



APPENDIX A: ELEMENTS OF C ab AND d ab 



Here we show the elements of C ab and d ab explicitly 

2(l + |ap) (1 + \a\ z Y 

(O12 = ^0^ (a*\u a ){v a \a*)(a*\v b ){u b \a*) + ^ ~f + Z^T^ <\\u a )(v a \a*){a*\v b ){u b \V, 

(c ab hi = ^^ai«-x«aU>ai«6><«6it), 

(O22 = fe^S (TI«a)K|a*)(a*li; t )(u 6 lt), (Al) 
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1 - e~ 7T M 2 (l — e ~ 7T — 7Te~ 7T ) 
^ Q ' > ' )l = 2(l + \a\ 2 ) ( ' a *^ Ua ^ Va ^ P ^ Vb ^ Ub ^ + (l + |of) 2 (TK>KIP(°)k>( u fclT>> 



^ = ferxrw(tl^)K|p(o)|^)Kit). (A2) 

2(1 + \a\ 2 ) 



APPENDIX B: DERIVATION OF THE RELATION (T50)) 

It is not difficult to perform the first, say, two summations over kg and fc^_i in the lhs of ([5 



k-l k-l-ki k-l-ki 

e e e 

fe,=0 fe 2 =0 fe;=0 



fc-^-fei fc«_ 2 / fc-(£-i)-fci fe € _! fc-(^-i)-fei fe<_r 

fel I ^ I ~t±l 



x l ■ ' ■ x t—l 



xi - xe+i xg+i - xi 



E- E 

fei=0 fci-i=0 

k-i-ki fcf-3 / fe-(^-2)-fci fe f _2 k-{i-2)-k 1 k e -2 



= E- E 

fc 1= k» 2 =0 



(xg-i - xe+i)(xt-i ~ xi) {x t - xi+i)(x e - xt-\) 

k-(i-2)-k! k e _ 2 



(Xi+i ~ X t )(x e+1 - Xg-i) 



(Bl) 



These expressions are already quite suggestive. Actually we easily perform the summation of the following form to 
get 

k—£ — ki fem— l / r-k—m—ki k m 



E 



„k m I X m+1 



m \(x m+ i - X £+ i)(x m+1 — Xi)-- - (x m+ i - X m+2 ) 

k—m — ki k ri 

X l+1 



{x£ + i - xe)(xe+i - Xi-i) ■ ■ ■ {xi+i ~ x m +i) 



-X 



( r e ~ m 



(Xm+1 - Xt+l) ■ ■ ■ {Xm+1 - X m ) {xt+l - Xl) ■ ■ ■ (xe+1 - X m ) 
k—(m—l)—ki k m -i k— (m— 1)— fci km— 



(x m +i - xi+i) ■ ■ ■ (x m +i - x m ) {xe+i - Xi)-- - {xg+i - x m ) 

The quantity in the square brackets in the first line of the rhs can be written as 



(B2) 



Jt-m 



(x m +i - x e+ i) ■ ■ ■ (x m +i - x m ) {xt+i - Xi)- - ■ (xi+i - x m ) 

1 



{Xi Xj ; 

e+l>i>j>m 



n ( x i- x i)-4- m j\{x i - Xj )+---+{-iy- m x t -? 1 n ( Xi - Xj ) 

i>j i>j i>j 



(B3) 

Since there can be no singularities at Xi = Xj ^ m) and the quantity in the last square brackets is a polynomial 
of (£ — m)th order all in xt+i, . . . , x m +i, one can easily deduce its form as 

f(x m )(xt+i - xi)(x l+ i - xt-i) ■ ■ ■ {xi+i - x m +i) 

x (x e - xt-i)(xe - xe- 2 ) ■■■{xi- x m+ i) 

). (B4) 

The remaining function f(x m ) can be fixed, for example, by the behavior around x^+i = x m , to be f(x m ) = 
(— l) e ~ m x m ~ rn . This means that the first term in (|B2[) is also expressed as 

p( \ k— (m— 1)— fci fcm— l 

- X k ~ l ~ kl fc m-l+l I\ x m) _ Xm 

II (Xi — x m ) [x m — X( + i)(x m — Xi) ■ ■ ■ (x m — x m+ i) 

i+l>i>m 
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The validity of (150)) is now evident. 

APPENDIX C: DOMINANT TERMS IN THE SUMMATIONS ([M]) 

The fact that the dominant contributions in (|54j) are due, in general, to those terms with the same set of indices 
{ij} = Wibi} = {0,2^2} = • • • = {a^+i&^+i} is understood by observing that the expression (|54p is nothing but the 
result of the summations of the terms with different powers of complex numbers and that no cancellations are possible 
among those terms with no relative phase, while strong cancellation are expected to occur among other terms with 
different phases. Indeed, (|54[) is nothing but the result of the summations of the form 

l—kki k e k—i—ki k e 

k—t—k\ ke- 1 

x N - 1 - e (x 1 /x e+1 ) k ^ ■ ■ ■ (x e /x e+1 ) k ' (x e+1 /x) k - e 

k e =0 

■ ■ ]T x N - 1 - e (x 1 /xf(x 2 /x 1 ) k ^ ■ ■ ■ (x l+1 /x e ) k K (CI) 

the absolute value of which is expected to be maximum when all complex variables x, Xi, . . . , X£+i are in phase. Since 
in the present case the absolute values of A an b n s are bounded |A a?i bJ < |An| = |Ai| 2 (we assume that there is no 
degeneracy |Ai| > | A2 1 and that the eigenvalues Xi have a nonvanishing relative phase), dominant contributions come 
from those terms with x\ = ■ ■ ■ = xg+i = An when x = An or x = A22 and otherwise with x\ = ■ ■ ■ = x^ + \ = x. We 
can thus easily extract the dominant contributions. Notice that when x = x\ = ■ ■ ■ = xi+\, the above summations 

N-l k-t k-l-kx fcf-i N-l-l N-l-l-k N-l-l-k-k t fc € _! 

EE- E = E E ••• E E-E ( C2 ) 

k=l fei=0 fcf=0 fc=0 fci=0 k e =0 a<k,k 1 ,...,k, +1 <N-i-e 

k + k 1 -\ \-k e+1 = N 

express nothing but the number of distributions of TV — 1 — I identical balls into I + 2 boxes allowing empty boxes, 
that is, jvCV+i. Furthermore, k = £ case is exceptional for in this case the summations over all and bj in (|49|) are 
trivially done. 



JV-l k-l fe-f-fei k—t—k-i k t 



EE E 

k=l fci=0 fe 2 =0 



N- 



E 

k e =0 
N-l k-l k-l-ki 

EE E 

k=l fci=0 k 2 =0 
N-l-e k' k[ 

E EE' 

k'=0 k'=0k'=0 
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